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Abstract. Some results on the ordered statistics of eigenvalues for one-dimensional 
random Schrodinger Hamiltonians are reviewed. In the case of supersymmetric 
quantum mechanics with disorder, the existence of low energy delocalized states 
induces eigenvalue correlations and makes the ordered statistics problem nontrivial. 
The resulting distributions are used to analyze the problem of classical diffusion in a 
random force field (Sinai problem) in the presence of weakly concentrated absorbers. 
It is shown that the slowly decaying averaged return probability of the Sinai problem, 
P(x, t\x, 0) <~ ln~ 2 t, is converted into a power law decay, P(x, t\x, 0) <~ t~^ 2p ^ 9 , where 
g is the strength of the random force field and p the density of absorbers. 
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1. Introduction 

An ordered (or extreme value) statistics problem can be denned as follows : given 
a set of Af ranked random variables X\ < x 2 < ■ ■ ■ < x n < ■ ■ ■ < xjy, one asks for 
the probability density W n ^(x) for finding x n at x. When the random variables are 
independent and identically distributed (i.i.d) according to a given distribution p(x), 
the distributions W n ^{x) may be easily related to p(x). The interesting point emerges 
when considering the Af — > oo limit : in this case, up to a rescalling x n = a n ^j + b n ^ y n , 
that depends on the original distribution p(x), the distribution of the rescaled variable 
y n converges to an universall&w w n (y). Three universality classes exist, corresponding 
to the nature of the tail of p(x), as demonstrated in the pioneering works of Frechet 
[7J and Gumbel (see the textbook [ID])- When p(x) has a tail of exponential type, one 
obtains the famous generalized Gumbel laws w n (y) = ex P n {y ~ eV )- The study 

of extreme events is relevant in many areas like meteorology, finance, computer science, 
statistical physics, etc, and in general the random variables of interest may not fulfill the 
hypothesis of statistical independence (see for example Ref. [18] or [1]). In this case the 
problem is much more difficult to analyze and no general results are known such as for 
i.i.d random variables. This situation occurs when considering the set of eigenvalues of a 
random operator, since eigenvalues are a priori correlated. For example, the distribution 
of the smallest eigenvalue of large random matrices has been obtained in a famous work 
by Tracy and Widom [271 [28] for the usual Gaussian ensembles ; other ensembles have 
been recently considered in Refs. [T31 HU 1221 ETJ- In the present paper, we will consider 
the case where random variables are eigenvalues of a one-dimensional (ID) Schrodinger 
operator with a random potential. 

The paper is organised as follows : in section [2] we define the problem and 
review some results. In particular we will focus ourselves on the case of disordered 
supersymmetric quantum mechanics for which we obtain a set of non trivial distributions 
w n (y). Applications of these results are discussed in sections [3] and [i] : in the context of 
ID disordered quantum mechanics and in the context of classical diffusion in a random 
force field with absorption. 

2. Ordered spectral statistics 

Let us consider the Schrodinger operator i? sca iar = — ^ + V(x) acting on functions 
defined on the interval [0, L], satisfying Dirichlet boundary conditions </?(0) = <f(L) = 0. 
We denote by E\ < E 2 < ■ ■ ■ < E n < ■ ■ ■ the (infinite) set or eigenvalues. We will 
be interested in situations where the potential V(x) is random. The section is devoted 
to the analysis of the probability density for finding the n-th eigenvalue E n at level 
E, denoted W n (E; L) = 5(E — E n ), where 7TT denotes averaging with respect with the 
disordered potential V(x). The knowledge of this set of distributions provides some 
spectral informations much more precise than the averaged density of states (DoS) 
g(E; L) = J2^=i W n (E; L), measuring the probability to find any eigenvalue at level E. 
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The determination of the distributions W n (E; L) has been achieved for different types 
of random potentials, what we review in the rest of the section. We will consider the 
limit L — > oo (in practice L must be larger than a certain scale charcateristic of the 
disorder) analogous to the J\f — > oo limit of the introduction. 

2.1. Fully localized spectrum 

The distributions W n (E; L) were determined when V(x) is a Gaussian white noise : 
the distribution of the ground state energy (n = 1) was obtained in Ref. [19], a 
result generalized to all eigenvalues in Ref. [23J. The case of a potential describing 
the superposition of repulsive scatterers at random positions, V(x) = J2i a i^( x ~ x i)i 
was considered in Ref. pJ] ; are uncorrelated and uniformly distributed on [0, L] 
with a mean density p ; the positive weights cti can be chosen random or not. These 
two models of random potential present different spectral densities and localization 
properties. They share however two facts [16J : (i) a suppression of the low energy DoS 
(Lifshits tails), compared to the free DoS, and (ii) a decreasing localization length as E 
decreases (low energy states get more localized). This explains the origin of the result 
obtained in Refs. |23] : after an appropriate model dependent rescaling of energies, 
y n = (E n — E n )/y/vax(E n ), the distributions W n (E; L) coincide, in the limit L — >• oo, 
with the generalized Gumbel laws w n (y) aforementionned for uncorrelated variables. 
The absence of energy level correlations is explained by the strong localization of low 
energy states trapped in deep potential wells [20J. 

2.2. Spectrum unfolding 

An interesting observation is related to the structure obtained in Ref. [23] for local- 
ized eigenstates in the L ^ oo limit : W n {E;L) ~ LN\E)j^y [LN(E)} n ~ l exp —LN(E), 

where N(E) = lim^oo \ J_ dE' g(E'; L) is the integrated DoS per unit length 
(IDoSpul) for L = oo. The change of variable £ n = LN(E n ) corresponds to unfolding 
the spectrum : it relates a set of random variables {E n } distributed with a non-uniform 
density g(E; L) ~ LN'(E), to random variables {£ n } distributed with a uniform density 
equal to unity. The fact that these latter variables are distributed according to a Poisson 
law ro n (£) = ~^zjy^ n ~ 1 demonstrates the absence of energy level correlations [20] . 
since m n (£) is also the probability that n — 1 eigenvalues lie in the interval [0, £]. 

2.3. Supersymmetric case 

The results reviewed in the previous subsection apply to the generic situation where 
all eigenstates are localised. However this does not exhaust all possible scenarii, as there 
exist one-dimensional disordered models exhibiting delocalized states. Such an example 
is provided by the supersymmetric quantum Hamiltonian 

#susy = + <P(X) 2 + <fj{x) = Q ] Q (1) 
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where Q = — ^- + 4>{x). This Hamiltonian appears in many interesting contexts (see 
Refs. [3 HI [25] for reviews). When 4>(x) has short range correlations and <p(x) = 0, we 
obtain low energy properties opposite to the one mentioned previously for -ff sca iar : (i) 
the disorder increases the low energy DoS (Dyson singularity) and (ii) the states get 
delocalized near the spectrum boundary at E = (this was demonstrated in Ref. [2] 
when 4>(x) is a Gaussian white noise). The derealization of the first eigenstates are 
responsible for energy level correlations and the distributions W n (E; L) are no longer 
expected to coincide with the generalized Gumbel laws. We now explain the principle of 
the method, introduced in Ref. [23J , allowing for the determination of these distributions, 
when 4>(x) is a Gaussian white noise such that <fi{x) = and 4>(x)(f)(x') = gS(x — x'). 

The starting point consists in converting the Sturm-Liouville (spectral) problem, 
H susy (p(x) = Eip(x) with boundary conditions (p(0) = <p(L) = 0, into a Cauchy problem, 
H S usyip(x', E) = Ei/j(x] E) with boundary conditions ip(0; E) = and ^'(O; E) = 1. The 
former admits solutions only for energy in a discrete set E e {E n }, while the latter has 
solutions WE G KL We introduce the notation £i, £ 2 , etc, for the lengths between the 
consecutive nodes of the wave function (Fig. [T]), and P(£) their common distribution, 
that will be obtained by studying the statistical properties of ip(x;E), or rather of 
some related Riccati variable. Then the probability for the eigenenergy E n to be at 
level E coincides with the probability for the n-th zero to coincide with the boundary 
Em=i £ rn = L, i.e. W n (E; L) oc (P * ■ ■ ■ * P){L), as illustrated on Fig. [l| This is the 
essence of the node counting method, also denoted phase formalism in the context of 
disordered systems [16J (or Dyson-Schmidt method for discrete models [T7]). 




4 ia 

Figure 1. Solution of the Cauchy problem for different energies, below, equal and 
above the ground state energy. 



The building brick is now the distribution P(£) and we explain how it can be 
calculated. Introducing the Riccati variable z = ifi'/if) — <f), we map the Schrodinger 
equation onto a Langevin like equation 

4-z(x) = ~E- z(xf - 2 z(x) (j)(x) . (2) 
dx 

In this language the distance £ n separates two consecutive divergencies of the Riccati 
variable, e.g. z(0) = +oo and z(£i) = — oo. The study of P(£) is therefore mapped onto 
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a first passage problem for a stochastic Markovian process, a problem solved by well- 
known techniques. We introduce the auxiliary "time" C z needed by the process z(x) in 
order to reach — oo for the first time, given that it has started at z. The characteristic 
function h(z; a) = e~ aCz obeys a backward Fokker-Planck equation 

B z h(z; a) — ah(z; a) (3) 

for boundary conditions h(— oo; a) = 1 and h'(+oo; a) = 0, where B z = 2gzj^z-^ — (E + 
z 2 )-^ is the generator of the diffusion. Finally the distribution P(£) can be obtained by 
an inverse Laplace transform 

POO 

h(+oo;a) = / d£P(£)e- a£ . (4) 



o 



An approximation scheme was proposed in Ref. [23] in order to solve Eq. ^ in the 
low energy regime E <C g 2 . This supposes that the distributions W n (E; L) are mostly 
concentrated on such energy scales, what is expected to occur when L — y oo. As a 
result, it was shown that 

/i(+oo;a)~ — - 1 = (5) 

cosh 2 y/a/N swsy {E) 

where 1/ £ = N smy (E) ~ 2g/ln 2 (g 2 /E) is the IDoSpul of the supersymmetric 
Hamiltonian ([I]) in bulk (L = oo). After spectrum unfolding £ n = LN susy (E n ), the 
distributions of the eigenvalues W n (E; L) = LN!, nsy (E) zu n (LN susy (E)) are given by 

dq e q S 
2m cosh 2 ™ y/q 

where is a Bromwich contour. The fact that these distributions deviate from the 
Poisson distribution (i.e. W n (E; L) differ from the generalized Gumbel laws) signals 
spectral correlations induced by derealization. These distributions were explicitely 
calculated in Ref. [25] by a generating function method : 



-«(£) = / S 77^= (6) 



The large £ behaviour is m n {£) ~ (2n-i)! ^ 2 "~ 1 ex P ~T"£" ^ e ^ rs ^ distributions are 
represented on Fig. [2j 

The sum of the distributions characterizes the averaged DoS 

oo 

mI) = J2Wn(E;L) (8) 

n=l 



for a finite interval. It coincides to the bulk density of states g(E; L) ~ LN' snsy (E) 
when eigenstates are localized on scale smaller than L and do not feel the boundaries. 
A log-normal depletion, g(E;L) ~ eI/jl ex P ~2gZ ^(d 2 /^)^ * s however obtained at 
low energy E <C g 2 exp —y / 2gL, illustrating the sensitivity of low energy delocalized 
states to the boundaries. The DoS is represented on Fig. [2] (dashed line) after spectrum 



unfolding. Conversely the boundary sensitivity of the averaged DoS g(E; L) may be 
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used as a localization criterium, what gives the energy dependence of the localization 
length f E ~ MnV/£)- 

Details on the results reviewed in this section can be found in Refs. [23j EE] • 



e * (localised case) 

IT r 








CO 

n =1 













Figure 2. Left : The distribution may &e interpreted as the distribution of the 

rescaled distance i between consecutive nodes of the wave function P(£) — vj\(l/t)/t ; 
the low energy wave function of the SUSY Hamiltonian therefore exhibits repulsion 
between nodes. For reference, we have also plotted (dotted black) the exponential 
distribution (absence of repulsion between nodes) characterizing localised eigenstates, 
the case discussed m § 2.1 & 2.2. Right : Distributions m n (£) of the five first rescaled 
eigenvalues = LN susy (E n ) . The sum (dashed line) corresponds to the averaged DoS 
in a finite interval, after spectrum unfolding with respect to bulk DoS ; depletion near 
£ = is due to the derealization. 



3. Supersymmetry broken by disorder 

As we have mentioned, the random Hamiltonians H sca \ ar and H snsy present opposite 
spectral and localization properties. This observation has led to question the interplay 
between the two types of disorder and consider the mixed case 

d 2 

H = -— + <j)(x) 2 + <j)'(x) + V(x) . (9) 

The case where 4> and V are two Gaussian white noises was analyzed in Ref . [UJ , however 
it turns out that the case where the potential V describes a random superposition of 
repulsive scatterers at random positions, V(x) = J2i a i$( x 

be related to 

an interesting problem in the context of classical diffusion (discussed in next section). 
In this latter case, the spectral density of the Hamiltonian ^ may be obtained from 
the distributions W n (E;L) obtained in the previous section by using a Lifshits like 
argument [23]. Let us denote by N(E) the IDoSpul of the Hamiltonian (|9]). In the 
limit of large weights — > oo (in practice ^> p & g) the impurities decouple the 
intervals free of impurity and impose Dirichlet boundary conditions at their locations. 
Thus the spectrum of ^ is given the addition of spectra of H susy on all intervals and 
we obtain the DoSpul N'(E) for the mixed disorders by averaging the DoS @ of the 
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supersymmetric Hamiltonian for a finite length 



N'(E)~p(g(E;L)) L , (10) 

where (• • -)l denotes averaging with respect to the length of interval, namely with an 
exponential law pe~ pL characterizing uncorrelated positions Xi. Using (6j8 ) we get 

N{E) ~ ; 9 . (11) 

sinh 2 ^p/N snsy (E) 

Above the threshold energy E c = g 2 exp — ^2g/p, the IDoS coincides with the one of the 
supersymmetric Hamiltonian. Below the threshold, the Dyson singularity is transformed 
into a power-law singularity 



A 4 "(fj ' (12) 

This power law behaviour has been very well confirmed by numerical calculations [24]. 
Moreover the numerics has shown that the exponent ^2p/g seems more general than 
the range of application of the Lifshits argument suggests, what we will comment on in 
the last section. 



4. Sinai diffusion with weakly concentrated absorbers 

The study of classical diffusion in a random environment has attracted a lot of 
attention in various contexts ranging from mathematics, statistical physics (glassy 
dynamics or polymer physics) and even finance. It may be described within the 
framework of Fokker-Planck equation d t P(x;t) = d x (d x — 2<p(x))P(x;t), where <ft(x) 
is chosen to be a Gaussian white noise with <f>(x) = pg and (f)(x)<p(x') — <p(x) 4>(x') = 
g8(x — x'). The model exhibits a rich dynamics as a function of the drift p [2J. In the 
absence of drift, for p = 0, the diffusion is controlled by overcoming of potential barriers 
and is extremely slow x(t) ~ ^\n 2 (g 2 t) (Sinai diffusion |22j). 

The question I would like to discuss now is : how does the introduction of absorbing 
sites affect the Sinai dynamics ? Let us go back for a moment to the free diffusion : in 
this case the return probability decays in time as a power law P(x, t\x, 0) = 1/ y/Ant. The 
dynamics is slowed down by the introduction of a random force field, what is reflected 
in the return probability by an extremely slow decay P(x,t\x,0) ~ g\n~ 2 {g 2 t) 1 related 
to the behaviour aforementioned. On the other hand, when a weak concentration p of 
efficient absorbers with large absorbing rates a ^> p is introduced, the free power law 
decay is transformed into an exponential decay P(x,t\x,0) ~ exp — 3(^p 2 t) 1//3 (Lifshits 
tail). The question is : what is the behaviour of the return probability when both 
random force field and absorbers are present ? Is the return probability increased or 
decreased ? 

The answer to this question can be obtained by a simple mapping to the model 
analyzed in the previous section. The effect of absorbers on the Sinai diffusion may be 
accounted for by adding to the Fokker-Planck equation a term : 

d t P(x; t) = d x (d x - 2<j)(x))P(x; t) - V(x) P(x; t) (13) 
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where V(x) = J2i a i $(. x ~ x i) describes absorbers located at a set of positions 
{xi}. The positive coefficient ctj measures the efficiency of the absorber located at 



Xi. The Fokker-Planck equation (13) may be mapped onto the Schrodinger equation 
—dt^(x;t) = H^(x;t) for the Hamiltonian ^ thanks to the transformation P(x;t) = 
ip{x;t) exp J x dx' 4>{x'). It follows that the return probability can be related to the 
spectral density (DoSpul) p(E) = N'(E) of Hamiltonian (|9| by 



P(x,t\x,0) = I dEp(E)e- Et . (14) 



We immediately deduce the large time behaviour of the return probability from (12) : 



P{x,t\x,0) -rv^for t»t c , (15) 



where the scale t c = 1/E C = exp ^2g/p is the time needed by the particle moving in 
the random force field to reach the closest absorber x(t c ) ~ 1/p where x(t) ~ - \n 2 (g 2 t). 
The power-law can be mostly explained by analyzing the lowest absorption rate 
E\[(j)(x),L\ in a finite interval : the fact that its distribution presents a log-normal 
behaviour Wi(E; L) ~ i exp — ^ ln 2 (g 2 /E) for E — > explains the exponent. The 



behaviour (15) was first obtained in Ref. [21] and later confirmed in Ref. [15J by the real 
space renormalization group method. The return probability decays slower than in the 
free case for a sufficiently weak concentration of absorbers p < g/8, whereas it decays 
faster for p > g/8. Note however that the origin of the decay in the present model 
is mostly due to absorption, as demonstrated in Ref. [15J by analyzing the survival 
probability which was shown to present the same power law decay up to a logarithmic 
correction S Xo (t) = J dx Pjx, t\x , 0) ~ (hit) rv 7 ^. 



5. Concluding remarks 

In this brief review, I have analyzed the ordered statistics (extreme value statistics) 
for eigenvalues of a random supersymmetric Hamiltonian. The fact that this quantum 
Hamiltonian possesses delocalized low energy eigenstates makes the problem nontrivial 
and escape from the universal generalized Gumbel laws obtained for generic ID random 
Hamiltonians with a full spectrum of localized states. 

These results have found some application in order to analyze the effect of dilute 
absorbers on the classical diffusion in a random environment (Sinai diffusion). We 
have obtained that the return probability presents a power law decay, mostly due to 
absorption. This conclusion relies on the analysis of the spectral density for a mixed 
random Hamiltonian with a "supersymmetric part" <fi 2 + <ft', where <f>(x) is a Gaussian 
white noise, and a "scalar part" V(x) = J2i a i <K X — x i)- Using a Lifshits like argument, 
it was shown that the DoS of this Hamiltonian presents a power law singularity at low 
energy p(E) = N'(E) ~ ^-i+v 7 ^?. The result is interesting by itself since only few 
power law spectral singularities are known in the context of ID Anderson localization : 
(i) the Halperin singularity, for randomly dropped attractive impurities of fixed weights 
Oi j = -v P| [IB] : p(E) ~ \E + !?;2|-i+2p/« for E ^ The pQwer law 
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behaviour for the supersymmetric Hamiltonian H svtsy for a finite 4>(x) = pg 7^ for 
which p(E) ~ E~ 1+fl [2]. 

Surprisingly, numerical simulations of the model ^ revealed that the exponent 
appearing in the IDoSpul (12) is much more robust than one would have expected on 
the basis of the heuristic Lifshits argument. Examining this observation has led us to 
a solvable version of the model by considering random weights on. Let us consider the 



Hamiltonian (|9j) for a finite drift <j)(x) = pg 7^ and with exponentially distributed 
positive weights = - exp — ^. This model becomes solvable for a particular value 
of the parameter v = ~ai : iiv = 2g the effect of the "scalar term" V(x) can be absorbed 
in a redefinition of the parameter p 



A*— n/=^* a + y (16) 

The IDoS of @ can be obtained by performing this substitution in the IDoS of H susy 
(that can be found in Ref. [2]), therefore : 

N{E) = % — 1 — , (17) 

^J v {VE/gf + N v {^E/gY' 

where J v {x) and N v {x) are Bessel functions [8]. Using their asymptotic behaviours we 

obtain the low energy IDoS N(E) ~ fly? (4^)^ ^ or ^ ^ 9 2 - Therefore we have recovered 

the power law with exponent Jlpj g (setting p — 0) within a solvable model, confirming 

the robustness of this exponent valid for arbitrary value of the ratio p/g. It is worth 

emphasizing that the solvable model applies to a regime with weights ol% ~ g while the 

analysis of the previous sections rather relied on the hypothesis that ojj 3> g, p. This 

result also generalizes the analysis to finite drift p 7^ (note that the same exponent 

\/ p 2 + 2p/g was obtained by Le Doussal in Ref. [15J by another method, in a different 

regime aj — > 00). Details of this analysis will be published in a forthcoming article [26J. 
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Appendix A. Diffusion constant of a single rf-dimensional Brownian curve 



In Ref. [25], the distribution of the maximal height of a Brownian excursion was 
shown to be related to the distributions (6[7), or more precisely to their sum ^^ =1 ^ n (0- 
In this appendix, we point out that the distributions 

J m 2m cosh y/q 
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generalizing m n (£) = ff2n(0 and considered in Ref. [25J, are relevant in order to 
characterize a certain property of the d- dimensional Brownian motion. 

Let us consider a <i-dimensional Brownian curve r(r) with r G [0, t] starting from 
the origin r(0) = 0. This curve is weighted according to the Wiener measure 

Vf{r)e-TD^ AT ^) , (A.2) 

where D is the diffusion constant, usually defined by D = lim^oo ^E,(f(t) 2 ), where 
E (■ ■ ■) denotes the expectation with respect to the Wiener measure ( A.2[ ). We now ask 
the question : is it possible to obtain an estimation of the diffusion constant from a given 
realization of the Brownian path by a time average ? For this purpose we introduce 
the "trajectory-dependent diffusion constant" 3i t = ^ J ^-^*(t) 2 , or preferably the 
dimensionless Brownian functional 

*mi-§^/t ? m 2 - < A - 3 » 

Note that this problem was also considered in Ref. [3] , starting from a different definition 
of the trajectory-dependent diffusion constant J ^r(r) 2 . 

The characteristic function of the functional x[r] may be written thanks to a path 
integral 



T{p) d = f E(e- p *^]) = [dR f rt) Vf(T)e-^^Hrr+^(^} {AA) 

J Jr(0)=0 



We recognize the integral of the imaginary time propagator (R \e~ tH \ 0) for a quantum 
harmonic oscillator with mass m — > 1/{2D) and pulsation uj — > \\[\- Using its well- 
known expression jH] we straightforwardly obtain 

r(p) = — d ) ^ ■ (A.5) 

cosh 2 l2y/p/d 



We check that E(x[r]) = 1 as it should. The variance is given by Var(x[r]) = ^ and 
vanishes in the limit of large dimension, as expected from the central limit theorem. 
In the case of even dimensions, we recognize the Laplace transform of the distribution 
(A.l). Therefore we deduce that the distribution of the functional x[^\ i s 

Pdix) = - 1 tX ) f° r dimension d even . (A. 6) 



4 2^4 

Expressions of the distributions 7r n (£) can be found in Ref. [25 
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